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ABSTRACT 

Strange-looking dust cloud around asteroid (596) Scheila was discovered on 
2010 December 11.44-11.47. Unlike normal cometary tails, it consisted of three 
tails and faded within two months. We constructed a model to reproduce the 
morphology of the dust cloud based on the laboratory measurement of high ve- 
locity impacts and the dust dynamics. As the result, we succeeded in the re- 
production of peculiar dust cloud by an impact-driven ejecta plume consisting of 
an impact cone and downrange plume. Assuming an impact angle of 45° , our 
model suggests that a decameter-sized asteroid collided with (596) Scheila from 
the direction of {aim, ^im) = (60° , -40° ) in J2000 coordinates on 2010 December 
3. The maximum ejection velocity of the dust particles exceeded 100 m/s. Our 
results suggest that the surface of (596) Scheila consists of materials with low 
tensile strength. 
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Subject headings: comets: general - comets: individual ((596) Scheila) minor 
planets, asteroids — general 



Introduction 



In this paper, we suggest one plausible explanation for the peculiar dust cloud of (596) 
Scheila. 

(596) Scheila is a large asteroid (113-120 km in diameter) orbiting the Sun in the 



outer region of t he main belt with the orbital period of 5.01 years f lTedesco fc DesertI 12002 



Usui et al.ll201ll ). The orbital elements of the object are typical for the outer main-belt, that 
is, the semi-major axis, the eccentricity and the inclination are 2.93AU, 0.164, and 14.7°. 
An unexpected dust cloud of (5 96) Scheila wa s discovered on 2010 December 11.4 with the 
0.68-m f/1.8 Schmidt telescope (ILarson Il2010l ). The observation with the same instruments 
on 2010 December 3.4 showed a diffuseness at r nagnitude 13. 2, about 1.3 mag bright er than 
that of the observation in the previous month (ILarson Il2010l ). Ilshiguro et al. I (120111 ) found 
that the dust particles ranging from 0.1-1 /xm to 100 /xm were ejected impulsively on 2010 
December 3.5 ±1.0 through the synchrone a nalysis of exte nded dust structure appeared 
after 2011 February. It is therefore likely that ILarson I (120101 ) observed (596) Scheila imme- 
diately after the dust emission. The total mass of the ejecta was estim ated to be (0.2-6) 



X 10^ kg, depending on the assurned pa. r ticle size and the mas s density (jjewitt et all 12011 



I 



Hsieh et al. 12011: Bodewits et al. 



never been detected (jJewitt et al 



201 



2011 



Hsieh et al. 



2011 



Ishiguro et al. 1120111). To date, gas ern ission has 



Bodewits et al.l 1201 ll ) . Accord- 



ingly, it is natural to think that the comet-like a ctivity was triggered by an impact . The 



impactor diameter w as estimated to be 20-50 m (jJewitt et al.l 1201 ll : iBodewits et al.l 12011 



Ishiguro et al. II2OIII ) 



Similar evidence for asteroid-asteroid impact was reported in another object. The dust 
cloud of P/2010 A2 (LINEAR) was discovered on 2010 January 6, showin g not only a comet - 
like extended dust cloud but also a mysterious X-shaped debris pattern (jJewitt et al.ll2010l ). 
It has been proposed that the dust ejecta was created by the impact of a small objects in 
2009 February or March (about 10 months before the discovery), although it cannot be ruled 
out that th e asteroid's rotational spin-up resulted i n a mass loss that formed a comet-like 
debris tail (ISnodgrass et al.ll2010l : IJewitt et al.ll2010l ). 



The morphology of (596) Scheila's dust cloud was also mysterious, in that it consisted 
of three prominent structures. Figure [1] (a) shows the image taken on 2010 December 12 
at Ishigakijima Astronomical Observatory with 1-m telescope and a 3ch (the g', Rc, and Ic- 
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bands) simultaneously imaging system. In the image, three c omponents appear: the northern 



ta il, southern t a il, an d westward tail ( llshiguro et al. Il201l[ ). Similar structures were found 



m 



Jewitt et al. 



(1201 if ). To date, no studies have addressed the physical implication of the 
mysterious morphologies of impact-triggered dust clouds. In this paper, we attempted to 
reconstruct the observed morphology of (596) Scheila's triple dust tails on the basis of the 
impact hypothesis. The image taken on 2010 December 12 at Ishigakijima Astronomical 
Observatory was compared to the model which considers the laboratory measurement of 
high velocity impact and dust dynamics, and derive the best-fit parameters. 



2. Model Description 



As a beginning we would like to clarify difference in shape between normal comets and 
(596) Scheila. The morphology of (596) Scheila cannot be explained by the sublimation 
of ice. Figure [T] (b) and (c) show the results of model simulations on 2010 December 12 
performed under the assumption of comet-like dust ejection. We applied continuous dust- 
ejection models using the parameters of 2 38P/Read as an an alog of a main-belt comet that 
was activated by the sublimation of ice flHsieh et al.l 120091) fFigure [1] (b)), and those of 
22P/Kopff as an analog of a Jupiter-family comet ( Ishiguro et al.l 120071 ) (Figure [11(c)). In 



these models, it is assumed that dust particles are ejected in cone-shape jets that are radially 
symmetric with respect to the Sun-object axis with a half-opening angle of 45° (238P/Read) 
and 60° (22P/Kopff). We considered continuous dust emission from two months prior to the 
observation on 2010 December 12. In these model images, the dust cloud smoothly extended 
in an almost anti-solar direction. There is only one dust tail in these images. The observed 
image differs from these simulation images in that it consisted of multiple tails. 

Secondly, we show simple impulsive emission models in Figure [1] (d)-(i). In these 
models, we assumed the isotropic dust emission of 1 /xm and lO/xm-particles with different 
velocities ejected on 2010 December 3. We adopted the ejection day based on the previous 



studi es; the dust emission should have occurred on 20 10 December 3.5±1.0 fjlshiguro et al 



20 111 ) but before 2010 December 3.4 (ILarson 1 12010| ). It may seem at a first glance that 
there should be dust particles ejected with high terminal velocity {vtmi ~ 100 m/s or higher) 
because the rim diameters of models with the terminal velocity of < 100 m/s look smaller 
than that of the observed image. In addition, small grains (<lyum) should exist because 
the dust cloud was deflected toward the anti-solar direction by solar radiation pressure. For 
comparison, assuming the mass density of 1670 kg/m^ (equivalent to t he m ass density of 



the Tagish Lake meteorite iHiroi fc Hasegawa I (120031 ): IZolensky et al. I fl2002l )). the escape 



velocity from (596) Scheila is 55 m/s. The estimated terminal velocity is two times faster 
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than the escape velocity. 

To reconstruct the observed morphology of the (596) Scheila dust cloud, we considered 
a new dust emission model, described below. We assumed that the dust particles were 
ejected in two different forms from an impact point, i.e., a conical impact ejecta curtain and 
a downrange plume, which are commonly observed in oblique impact experiments. Figure 
|2] shows an example of the laboratory oblique impact experiment. It was conducted with a 
two-stage hght gas gun at the Institute of Space and Astronautical Science (IS AS), Japan 
Aerospace Expl oration Agency (J AX A). A 3.2-mm steel sphere was accelerated with the 



gun in a sa bot (|Kawai et al.ll2010[) to 4.521 km/s, which is typical impact velocity in outer 
main-belt f lO'Brien et al.ll201l[ ). Soon after the impact, the dust particles were ejected as 
the luminous downrange plume with high speed in the horizontal direction with respect to 
the local surface. Later, the conical ejecta curtain was grown. We simplified these forms, 
as shown in Figure |2](e). The conical impact ejecta curtain is symmetrical with respect to 
a vector normal to the asteroid surface {acone, Scone) with a half-opening angle of 6. Dust 
particles are assumed to be ejected between 9-A9/2 and 9+A9/2 (see the shadowed area in 
Figure [2] (e)). We modeled the downrange plume as a stretched cone with a central axis 
ranging from (ald^„, 51du,n) to (a2rf^„,, (52d^„). 

The dust particles ejected with the ejection velocity of Vej would decelerate by the 
asteroid's gravity and reach the terminal velocity Vtmi- We applied the power law function 
of the ejection velocity of dust particles: 

Vej = VoO^ 

where Vq is the reference ejection velocity (m/s) of the particles radius a=lxlO~® (m), k is 
the power index of size dependence of the ejection velocity, and M^qq and R^qq are the mass 
and radius of (596) Scheila. We considered the energy conservation for the terminal velocity 
Vtml in the second equation of Eq. (1). 

An power-law size distribution with index q was used. The number of dust particles 
within a size range of a and a+da is given by: 

,w X , r Noa'^ da for a™„ < a < amax 

i\[ajda = s J. (2) 

to for a < amin, a > amax 

where a^ax and amin are the maximum and minimum particle sizes, respectively, and Nq 
represents the reference dust production rate. We fixed amm=0.1 /im and g=-3.5 (discuss 
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later), because particles much smaller than the wavelength are inefficient scatterers in th e 
optical wavelength. The size distribution exponent was also fixed to g=-3.5 (jDohnanyilll969l ) . 
The maximum size amax can be derived when the size-dependent velocity v^j becomes equal 
to the escape velocity from (596) Scheila. 

The trajectories of the particles were computed from the terminal velocity and the ratio 
of the force exerted by the solar radiation pressure and the solar gravity (/3). It can be 
expressed as f3=KQpr/ pa , where K = 5.7x1 0"'^ kg m~^, and Qp r is the radiation pressure 



coefficient averaged over the solar spectrum (IBurns et al. Ill979l ). We assumed Qpr=l. p 



denotes the mass density of dust particles. We supposed that the mass densities of the 
dust particles was 1670 kg/m^. The model images were reconstructed using the Monte Carlo 
approach for the p arameters above (llshiguro et al.ll2007l : Ilshiguroll2008l : ISarugaku et al.ll2007l : 
Hsieh et al.ll2009l ). We calculated the positions of dust particles at a given time by solving 
the Keplerian equation. We considered the impulsive dust emission on 2010 December 3 as 
stated above. The free parameters of our model are listed in Table 1. 



3. Results and Discussion 

Multiple simulations are carried out using various parameter sets, and the resulting 
model images are then visually compared to the data to find plausible model parameters. 
As the result, we obtained the best-fit values. Figure [3] shows some example results for 
the conical ejecta curtain. First, we noticed that the southern tail and westward tail (see 
figure [1]) could be reproduced by a conical ejecta curtain given a half-opening angle of 50° . 
Our observations were consistent with the case in which the central axis at {acone, ^cone) 
of the downrange was (90°, -15°) in the J2000 coordinate system. Similarly, we examined 
the dependence of the downrange plume {adwn, Sdwn) on the central axis. We notice the 
integral along the great circle joining from {aldwn, Sldwn) = (150°, -^40°) to {a2dwn, S2dwn) 
= (180°, +50°) matches the observed image. As 6 increased, the opening angle of the 
ejecta became broader; as Vq increased, the dust cloud extended more widely. Accordingly, 
these two variables, 6 and Vq, were well determined by comparison to the observed images. 
We estimated ^=15° (the downrange plume), ^=50° and d6=10° (the conical curtain). 
We derived Vq = 190 m/s (the downrange plume) and Vq = 80 m/s (the conical curtain), 
suggesting that the maximum speed was 340 m/s (the downrange plume) and 140 m/s (the 
conical curtain) for 0.1- pm particles. As k decreases, the near- nuclear dust cloud becomes 
brighter. The images we captured were consistent with the simulation images for a value of 
k ~ —1/4. By combining the conical ejecta curtain and the downrange plume models, we 
obtained the best-fit image (Figure H]). Note that we succeeded in the reproduction of the 
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triple tails by single impact model based on the idea that is commonly observed in oblique 
impact experiments. 

As we mentioned above we fixed two variables: q and amin- If there are a large amount of 
s mall particles in the dus t cloud, it should look bluer by Rayleigh scattering. As we noticed 



m llshiguro et al. I ( 1201 ll ) , no significant difference appears in the morphology observed in 
three different optical channels. The observational evidence implies that the diffuse cloud 
consisted of dust particles large enough to scatter optical light (i.e., 27ra/A > 1, where A 
denotes the optical wavelength). The size distribution exponent q=-3.5 was applied because 
it is the typical to impact fragment ( iDohnanyil Il969l ) . We performed the fitting with g=-4 
and g=-3, but could not obtain the plausible result. Therefore, the initial assumptions of q 
=-3.5 and amm=0.1/im seems to be reasonable. 

In the laboratory experiments, the central axis of the conical curtain is usually per- 
pendicular to the local surface. The downrange plume appears along the trajectory axis of 
the impactor. Therefore, we can derive the impactor's trajectory if the impact angle with 
respect to the lo cal surface is known. Th e most probable impact angle on arbitrary plane- 
tary body is 45°( lGault fc Wedekindlll978l ). If an impact angle of 45° is assumed, it is likely 
that a small asteroid collided with (596) Scheila from the direction of {aim-, ^im) = (60° , 
-40° ). This result suggests that the impactor collided with (596) Scheila from behind. The 
angle between the central axis of the conical curtain and that of the downrange plume in 
the best-fit model is 85°, which is potential value based on the impact experiments (Figure 
2(b) in this paper and Figure 19 of ISchultz et al.l (120071 )). The derived half-opening angle of 
50° is consistent with the results obtained for the rock y ejecta with velocity from hundreds to 
thousands m/s through the laboratory measurement ( iGault &: Heitowit Ill963l ). The power 
index of size dependence of the ejecta velocity {k ~ —1/4) is smaller than that of comets 
(i.e. k ~ 1/2, typica l of hydrody namical gas drag) but within the range of the laboratory 
impact experiments ( iGiblinI 119981 ). Our model predicts that up to 140/im particles could 
escape from (59 6) Scheila. In fact, 10 0/im particles were found in the observed images after 
2011 February (llshiguro et al. Il201ll ). The velocity of dust particles depends on the tensile 
strength of the surface materials when the impact process is dominated by the material 
strength rather than gravity. The maximum ejecta velocity for solid and porous targets 
measured in the laboratory exceed s 10 x (l^/p/j)"'^, wher e Yj; is the tensile strength and pt is 
the target density (in Figure 18 of iHousen fc Holsapple I ( 1201 ll ) and references therein), i.e., 
Vmax > 10 X {Yt/ PtT'^- Since the maximum ejecta velocity is 140 m/s for the conical ejecta, 
Yt < Pt{vm.ax/^Oy ~ 0.3 MPa. Our result on the ejecta velocity suggest that the surface on 
(596) Scheila was covered by materials with low tensile strength. 
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4. Summary 



So far we have outlined a plausible explanation for the peculiar dust cloud of (596) 
Scheila. We constructed a model of the morphology based on experiments of high velocity 
impacts. The values of the parameters were free and obtained from fitting the observed 
i mage on 2011 Decern ber 12. We found that the morphologies on 2011 December 17 and 19 
( llshiguro et al. Il201ll ) were also reproduced with the same model parameters. In summary, 
we find that: 



1. The morphology of (596) Scheila can be reproduced by an impact-driven ejecta plume 
consisting of an impact cone and downrange plume. 

2. The maximum ejection velocity of the dust particles exceeded 100 m/s. 

3. Assuming that an impact angle of 45°, the impactor collided with (596) Scheila from 
the direction of (60°, -40°) in J2000 coordinates. 

4. The surface of (596) Scheila consists of materials with low tensile strength (~0.3 MPa). 



With the previo us studies fjjewitt et al.ll201ll : iHsieh et al. Il201ll : iBodewits et al.l 12011 



Ishiguro et al. Il201ll ). we definitively conclude that a decameter-sized asteroid collided with 
(596) Scheila from behind on 2010 December 3. 



Research at Seoul National University was supported by the National Research Foun- 
dation of Korea and the Seoul National University Foundation Research Expense. This 
study was supported by ISAS/JAXA as a collaborative program with the Space Plasma 
Experiment. 
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Fig. 1. — (a) Image observed on 2010 December 12 at Ishigakijima Astronomical Observatory. 
(b)-(i) Images of cometary dust ejection models for the observation on 2010 December 12. 
(&)(h) Results of the continuous dust ejection rnodels using parameters from 238P/Read 
jHsieh et al.lboOQl ) and 22P/Kopff Jlshiguro et al.lboOTh . (d)-(i) The resuhs of an impulsive 
isotropic dust ejection. We assumed 1 /im-particles with the terminal velocities of (d) 57 
m/s, (e) 100 m/s, and (f) 200 m/s, and 10 /im-particles with the terminal velocities of (g) 
57 m/s, (h) 100 m/s, and (i) 200 m/s. In all panels, the emission source, (596) Scheila, is at 
the center of each image, and the field of view is 3.5' x 5.8'. 
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Fig. 2. — Sample images of the laboratory impact experiment. The experiment was con- 
ducted with a two-stage light gas gun at the Institute of Space and Astronautical Science 
(ISAS), Japan Aerospace Exploration Agency (J AXA). A 3.2-mm s teel sphere was accel- 
erated to 4.521 km/s with the gun using a sabot (IKawai et al.ll2010l ). The material of the 
target was serpentine. The chamber of the target was evacuated to a pressure of 7 Pa. The 
time resolution of the high-speed camera was 4 /is. The number at the bottom right of each 
frame shows the elapsed time, t, in microseconds, with t=0 at the impact, (a) An image 
before the impact, (b) Impact of the projectile forming a luminous downrange plume, (c) 
Dust particles were ejected from the target while the luminous downrange plume passed with 
high speed, (d) The growth of the conical ejecta curtain, (e) Schematic diagram showing 
how our simulation modeled the impact phenomena. Impact ejecta consist of a conical ejecta 
curtain and downrange plume. 
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Fig. 3. — Images of conical curtain ejection models for December 12, 2010 for 6=50° , 
A6'=10° , Vo=80 m/s, and different jet directions as labeled, where acone is constant for each 
row of models and Scone is constant for each column of models. In all panels, the emission 
source, (596) Scheila, is at the center of each image. 
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Fig. 4. — Result of numerical simulation reproducing the observed morphology on 2010 
December 12. In the model, we considered dust ejecta consisting of an impact cone and 
downrange plume. Filled circle indicates the asteroid's position; arrow indicates the anti- 
solar direction. The northern and southern tails were swept back by radiation pressure from 
the Sun, while the westward tail extended linearly because the orientation was close the 
anti-solar direction. 



- 15 - 



Table 1: Input and best-fit parameters for tlie dust ejection model 



Parameter Input values Best-fit values 

Conical curtain 

e [°] 10 - 90 with 5 interval 50 

56 [°] 1, 3, 5, 10, 15, 20, 30, 40 10 

ttcone [°] - 360 with 15 interval 90 

5^one [°] -90 - +90 with 10 interval -15 

Vq [m/s] 60 - 300 with 10 m/s interval 80 

k -1/6, -1/5, -1/4, -1/3, -1/2 -1/4 

q -3.5 (fixed) - 

amin [m] 1.0x10"'^ (fixed) 

amax[M defined as Wej(amaa;) =55 m/s 4x10"^ 
Downrange plume 

6 [°] 10 - 90 with 5 interval 15 

ttlrfion [°] - 350 with 15 interval 150 

5ldn,n [°] -90 - +90 with 10 interval +40 

(^'^dwn [°] 10 - 360 with 15 interval, a2dwn > cddwn 180 

(52rf^„ [°] -90 - +90 with 10 interval +50 

Vo [m/s] 60 -300 with 10 m/s interval 190 

k -1/6,-1/5,-1/4,-1/3,-1/2 -1/4 

q -3.5 (fixed) - 

amin [m] 1.0x10"''' (fixed) 

amax [m] defined as fej(0'max)=55 m/s 1.4x10"^ 



